Thermal transport in one-dimensional spin gap systems 
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We study thermal transport in one dimensional spin systems both in the presence and absence 
of impurities. In the absence of disorder, all these spin systems display a temperature dependent 
Drude peak in the thermal conductivity. In gapless systems, the low temperature Drude weight is 
proportional to temperature and to the central charge which characterizes the conformal field theory 
that describes the system at low energies. On the other hand, the low temperature Drude weight of 
spin gap systems shows an activated behavior modulated by a power law. For temperatures higher 
than the spin gap, one recovers the linear T behavior akin to gapless systems. For temperatures 
larger than the exchange coupling, the Drude weight decays as T~ 2 . We argue that this behavior 
is a generic feature of quasi one dimensional spin gap systems with a relativistic-like low energy 
dispersion. We also consider the effect of a magnetic field on the Drude weight with emphasis on the 
commensurate-incommensurate transition induced by it. We then study the effect of nonmagnetic 
impurities on the thermal conductivity of the dimerized XY chain and the spin-i two leg ladder. 
Impurities destroy the Drude peak and the thermal conductivity exhibits a purely activated behavior 
at low temperature, with an activation gap renormalized by disorder. The relevance of these results 
for experiments is briefly discussed. 



I. INTRODUCTION 

The past many years have seen a resurgence of inter- 
est, both theoretical and experimental in quasi-one di- 
mensional spin gap systems. Well known examples of 
systems with a gap are sain chains with dimerization, 
frustration and anisotropyacl. Another interesting exam- 
ple is the two leg spin S — \ ladder which was pro- 
posed as a toy model for the-pseudogap phase in high 
temperature superconductors. u Renewed interest in these 
systems was^ triggered by the availability of anisotropic 
materialsotm in which the magnetic properties of the 
insulating phase could be ascribed to one or quasi-one 
dimensional spin systems. The dynamical properties of 
these quasi-ld spin phases have been extensively stud- 
ied using standard techniques like neutron scattering and 
NMR measurements. OH More recently, heat transport is 
being used as a complementary probe to study low di- 
mensional spin systems. Measurements of thermal con- 
ductivity have been carried out in systems, pich as the 
spin chain materials SrCu02 and S^CuOsEtEj, the spin 
Peierls system CuGeQgiiaDd the spin ladder materials 
(Sr, Ca, La)i4Cu2404iE3ll3'O. The huge anisotropy seen 
in the thermal conductivity in the directions parallel and 
perpendicular to the chains or the ladders, indicates that 
magnetic excitations of these quasi- ID systems do play 
an important role in heat transport. This is further con- 
firmed by measurements in the presence of a magnetic 
fieldO. 



Various attempts have been made to extract from these 
measurements, the purely magnetic contribution to the 
thermal conductivity. This is often done by subtracting 
a phonon background calculated within a Debye model. 
However, in order to account for the entire magnetic con- 
tribution, one needs to understand the interactions of 
the spin excitations of the low dimensional spin systems 
with themselves and, with defects and phonons. This is 
a non-trivial problem since the spin excitations are not 
necessarily weakly interacting, and the form of interac- 
tion of these spin excitations with phonons or defects is 
usually rather complicated. Consequently, experimental 
results have been fitted using various phenomenological 
kinetic theory expressions for non-interacting spinons or 
magnons. This effort to obtain the purely magnetic con- 
tribution to the thermal conductivity has stimulated the- 
oretical studies of t hermal , transport in spin chain and 
spin ladder systems .li§li!30 In the absence of extrinsic 
scattering such as phonons or defects, some studiesE£l 
showed that the frequency dependent thermal conduc- 
tivity k(lo,T) = ttk(T)8(uj), where h(T) is the thermal 
Drude weight. However, the prude weight extracted 
from finite size zig-zag laddersO seems to be at odds 
with the idea of an infinite thermal conductivity in spin 
systems without disorder. 

In this paper, we use analytical methods to revisit the 
problem of the thermal conductivity for various quasi Id 
spin systems with special emphasis on the two leg spin-i 
ladder. In the absence of impurities, we present results 
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FIG. 1: Schematic representation of the thermal Drude 
weights in spin gap systems without impurities 



which should be valid for spin gap systems possessing 
low energy triplet excitations and gapless systems irre- 
spective of the details of the nature of the interaction. A 
schematic representation of the thermal Drude weight for 
a spin gap system is shown in Fig.0. We also study the 
effect of one/many impurities on the thermal conductiv- 
ity of the spin ladder and show that in the presence of 
impurities, the thermal conductivity is not simply given 
by the Drude weight times a temperature independent 
scattering time. 

The paper is organized as follows: in Sec. [n], we discuss 
the temperature dependence of the Drude weight for dif- 
ferent spin systems ranging from gapless integrable spin 
chains to spin gap systems like the spin ladder and the 
spin— 1 chain. In Sec. |l|, we discuss the effect of im- 
purities on the thermal conduction. In particular, we 
use the Landauer approachEfl to evaluate the effect of a 
single non-magnetic impurity on the thermal conductiv- 
ity of the ladder and the XY-chain. We then study the 
effect of a finite concentration of impurities on the lad- 
der. Finally, we present a comparison of our results to 
experiments and other theoretical work on the subject. 



II. TRANSLATIONALLY INVARIANT 
SYSTEMS 

In this section, we briefly outline the general definitions 
of the thermal current and the thermal conductivity cal- 
culated within linear response theory. We then use this 
formalism to calculate the dc thermal conductivity of var- 
ious spin systems with and without a gap to low energy 
excitations. The examples considered are: gapless inte- 
grable spin chains described by a conformal fixed point, 
the spin-i ladder, the dimerized XY chain and lastly the 
case of massive bosons. We also consider the effect of a 
magnetic field on the thermal Drude weight. 



A. Definition of thermal current and thermal 
conductivity 

We consider a system defined by a Hamiltonian density 
H(x) so that the total Hamiltonian is H = J dxH(x). 
Conservation of energy leads to the continuity equation 



where j e is the energy(thermal) current of the system. In 
the absence of charged excitations, the energy and ther- 
mal current are equivalent. Eq.(|l|) permits a definition 
of the thermal current in terms of Jthe Hamiltonian den- 
sity. Within linear response theoryuJ, the energy current 
response function at temperature T > reads: 



x{u,T) = dx 



dte % 



[? e (M),je(0,Q)]>, (2) 



where (. . .) indicates both quantum and thermal averag- 
ing. It is often easier to use the imaginary time formalism 
to calculate x(* w n) where uo n are the Matsubara frequen- 
cies and then analytically continue to real frequencies 
iuj n — > u> + iO to obtain x( w i T).l3 The frequency depen- 
dent thermal conductivity is then given bylij: 

k(cj,T) = ^-[ X (0,T)-x(u>,T)], (3) 

ILOL 

In general, in the absence of phonons or impurities, the 
total thermal current J e {t) = J dxj e (x,t) is conserved. 
This conservation permits an alternative but equivalent 
formulation of the thermal conductivity ([}]) 



K (w,T) = 



1 



2LT 2 



({J e (t),J e (0)}W wt dt, (4) 



where L is the system size. Since J e is conserved, 
total current is time independent and 



k(co,T) 



LT 2 



(J 2 e )6(co) = R(T)S(w) 



the 



(5) 



d t H(x,t) + d x j e (x,t) = 0, 



(1) 



This implies an infinite dc thermal conductivity, with a 
temperature dependent Drude weight H(T) which van- 
ishes at zero temperature. This thermal Drude weight 
has been studied numerically for some spin gap systems 
in Refs.[l5|,|l7]. In the following sections, we present an 
analytical discussion of the behavior of the Drude weight 
in various gapless and gapped quasi-one dimensional spin 
systems. 



B. Gapless Integrable Spin Chains 

In this section, we present results for the thermal 
Drude weight of integrable spin chains which are char- 
acterized by a vanishing singlet-triplet gap. One exam- 
ple of such a system, is the integrable spin-i Heisenberg 
model, which is known to have a Drude weight that vaH,- 
ishes linearly as temperature goes to zero k = 7r 2 uT/3cil, 
where v is the velocity of spin excitations. Other interest- 
ing systems, are the various integrable generalizations of 
the spipuji Heisenberg spin chain, like the snin— .S* chain 
modelsBB and SU{N) spin chain models^'EirEll This 
description allows one to easily obtain the low temper- 
ature thermal conductivity of these chains. The long 
wavelength behavior of these integrable systems are de- 
scribed by a conformal invariant fixed pointEZBEj. This 
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description allows one to easily obtain the low tempera- 
ture thermal Drude weight of these chains. The effective 
Hamiltonian of these systems has the generic form 



H = J dx[H R {x)+H L (o 



(6) 



where H. R and Tli describe right and left moving chi- 
ral modes. In addition, chirality imposes the constraints 
Ti. R (x,t) = 7i R (x — vt) and Hi(x, t) = Hi(x + vt). This 
leads to the following relation 

d t (H R (x,t)+H L (x,t)) = -vd x (H R (x,t) -H L (x,t))(7) 
which results in a thermal current density: 



J e = v / dx[: Hr{x) : - : Hl(x) :] 



(8) 



As before, since [H, J e ] = 0, the Drude weight is given 
by 



k(T) 



LT 2 



(9) 



Moreover, since there is no interaction between the right 
and the left moving modes, (J 2 ) = v 2 (H 2 ). One thus im- 
mediately obtains the result Kwz(T) — irC v (T)v 2 , where 
C v is the specific heat of these modes. For conformally 
invariant modes with a central charge c, the specific heat 



is given by C V (T) 
weight: 



7TT 

3v 



c, leading to a thermal Drude 



(10) 



For the integpaJsle spin-1 chain at the Takhtajan- 
Babujian pointE3£3, this weight can also be recovered 
from explicit calculations using the Majorana formalism 
to be discussed in the forthcoming sections. For a theory 
described by a free massless boson like the spin-^ Heisen- 
berg chain, which has a central charge c = 1, this weight 
is tt 2 Tv/3, which can also be checked by direct calcula- 
tions of the thermal susceptibilitidlZl For systems with a 
Luttingcr liquid like descriptioifj with a Luttinger ex- 
ponent K, the present derivation illustrates clearly that 
the weight R is independent of the Luttinger exponent 
or equivalently, the compactification radius of the free 
bosonic Luttinger field. Considering the case of XX Z 
chains, this result implies that the thermal Drude weight 
is independent of the anisotropy J z / J xy which is in agree- 
ment with Bethe Ansatz calculations, on the XXZ spin 
chain in the Luttinger liquid regime£il 



Spin- 



| ladder 



Here and in the following sections, we focus exclusively 
on spin gap systems. We first apply the formalism of 



Sec. U to the clean two leg spin ladder. The Hamiltonian 
of the two leg spin ladder is 



H = Jn 



J± ^ Sj,i • Si- 



(11) 



liP are spin-i operators, and the exchange 



where the S. 

constants Jn , J± > 0. For weak interchain coupling 
Jj_ <C J\\ , the spin ladder can be describecLby a contin- 
uum theory of spinless Majorana fermionsEU. The con- 
tinuum Hamiltonian reads: 



3 

H = / dxH a (x), 

n ** 



(12) 



n a (x) 



a=0 

-iv 

T 

<l{x)d x n{x 



[e R {x)d x e R (x) - m R )(x)c R (x) 

{d x e L )mi{x)]+irn a e R {x)a\$) 



where the velocity of the Majorana fermions v = ^ J\\a (a 
is the lattice spacing). Physically, the Majorana modes 
£ R i (a = 1,2,3) with masses 01^2,3 = J±/(2tt) = A 
describe triplet excitations with a gap A and £ R L with 
mass too = — 3J_l/(27t) = — 3A describe singlet excita- 
tions. We remark that the bosonized version of the low 
energy Hamiltonian ( |l2[ ) describes more geneijal gftin lad- 
der models than the one considered in (|ll])H3S Using 
M, the energy current for the ladder takes the form 



Je(x) = ^J e Q (x), 



(14) 



a=0 



ie R d x e R - (d x c R )c R + c L d x Ci - mim 



From ( |l4|) and (||), the total Drude weight for the spin 
ladder is found to be 



k(T) 



K a {T) — k (T) + 3k (T) 



(15) 



Since the Majorana fermions are essentially free, the cor- 
respondence between Majorana and Dirac fermions can 
be used to evaluate the Drude weight h a (T) 



k a (T) 



1 

8T 5 



dk- 



cosh' 



(16) 



and 



where the energy dispersion e a (k) — y (vk) 2 + 

A ~ — is the lattice induced ultra-violet cutoff. The 

a n 

details of the calculation are presented in Appendix W. 

The thermal Drude weight of the spin ladder is now given 
by 



k(T) 



3 
8T2 



4 1,2 



dk- 



v k 



-A 



cosh' 



y it 



1 

8T 2 



41,2 



dk- 



v k 



-A 



cosh' 



fl7) 



At low enough temperatures T <C J±/(2tt), the triplet 
excitations are the dominant carriers of heat and 



fi(T) 



ttA 2 
~2T 



■ ve 



— A/T 



(18) 
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For temperatures A <C T <C Jn, the coupling between 
the two spin half chains becomes irrelevant and we re- 
cover k oc T i.e., it is the sum of the P rude weights of 
two independent spin— i chains cf. Sec. II B. For tem- 



peratures T S> Jy, the temperature dependence in the 
integrand of ([It]) becomes negligible and since the k in- 
tegral is bounded on a lattice, the thermal Drude weight 
decays as k oc T~ 2 . The prefactor depends on the cut- 
off A and it is reasonable to assume that the continuum 
theory over-estimates this prefactor. To summarize, the 
Drude weight of the spin ladder has three regimes :i) 
at very low temperatures, T <C A, we obtain the expo- 
nential behavior ( |l8| ) ii) for intermediate temperatures, 
A < T < J,,, k oc T and hi) for T > J\\,R ~ l/T 2 . 



Since k{T — + 0) = 0, this implies the presence of at least 
one maximum in k at a hnite temperature for a lattice 
model. We expect k to have a peak in the vicinity of 
T ~ J|| (cf. Fig.|l|). We note that the numerical results 
for k(T) for the ladder presented in Ref.|l^ confirm our 
picture. 

To study the effect of an applied magnetic field h on the 
thermal conductivity, we first note that the effect of the 
magnetic field is to alter the dispersion of the triplet. The 
degenerate triplet dispersion ei(fc) now splits into three 
branches ei{k) + h,e\(k) and t\{k) — h and the singlet 
dispersion eo(fc) remains unaltered. The Drude weight in 
the presence of the field is now given by: 



J 



1 

8T 2 



dk 



foi(fc) 
dk 



cosh 



(ei(k) + h) 2 



r 



£l (fc) 2 



(siGh*) cosh 2 (Sl^+!t) cosh 2 (SM) 



1 

8T 2 



dk- 



; 4 k 2 



cosh 



2T 



(19) 



There are now two regimes of interest: h -C A and /t, > A. 
In the former case, the effective gap A — h dominates 
the thermal conductivity and k oc e" 
magnetic field leads to a sufficient enhancement of the 
low temperature thermal conductivity. The physical rea- 
son is that the increase of the number of triplet exci- 
tations with S z = +1 strongly dominates the diminu- 
tion of the number of excitations having S z = — 1. For 
h ~ A, thc dispersion of the Majorana fermions describ- 
ing the S z = +1 sector is no longer relativistic-like but 
quadratic, e(fc) oc fc 2 , resulting in k ~ T^/ 2 Finally, 
for h > A, the gap in spin ladder is closedL3L3, and the 
fermionic excitations have an effective linear dispersion, 
leading to k(T) = Z-^-v(h), where the effective Fermi 
velocity v(h) = vy/l - (A/h) 2 . 

Let us note that the above results (|l7|), (19) are also 
relevant for spin-1 chains. Indeed, spin-1 chains are also 
described at low energy by massive Majorana fermions, 
with a Hamiltonian similar to (|l^), except that the sin- 
glet mode £° is absentO. This mapping to massive Ma- 
jorana fermions originally derived for a spin-1 chain with 
bi-quadratic interactions in the vicinity of the Takhtajan- 
Babujian pointE3c3, is also expected to provide a qual- 
itatively description of the low energy properties of the 
Heisenberg spin-1 chain. Therefore, the thermal conduc- 
tivity of the spin-1 chain is easily obtained by taking the 
limit eo - > 00 in Eqs. © and (y_9|). For the spin-1 chain, 
the low temperature behavior of the Drude weight in the 
thermal conductivity is still given by (|l8| ) . The main dif- 
ference between the ladder and the spin-1 chain stems 
from the fact that while in thc former the gap to triplet 
excitations is small, in the latter the gap A is of the 
order of the Heisenberg exchange J (A = 0.41J). Con- 
sequently, the intermediate regime of linear temperature 



dependence of k can hardly be observed in the spin-1 
chain. However, reasonably strong bi-quadratic interac- 
tions can reduce the gap appreciably rendering an ob- 
servation of an intermediate linear regime possible. This 
predicted linear behavior might in fact be observable in 
the compound LiV2Ge06 which is expected to have size- 
able biquadratic interactionsEZI. 



D. Dimerized XY Chain 

We consider a spin-i XY chain with alternating ex- 
change in an external magnetic field h, described by the 
Hamiltonian: 

H = Jl(^2n^2n+l + ^2n^2n+l) 

n 

+ ^ JiiSinSzn-l + ^2n^2n-l) ~~ hS„ (20) 
n 

Using thc Jordan- Wigner transformatiorS, 

Sn = 4 cos I a m a m j , 
\m<n / 



nz f _J_ 

&n ~ a n a n „ • 



(21) 



where the a, are fermion annihilation and creation op- 
erators, the Hamiltonian (|2^) can be rewritten as: 

H = Ji 5^(4n+l°2n + H - C -) + J 2 ^( a 2n-l a 2n + H.c(J22) 

n n 

Diagonalizing the above Hamiltonian, we obtain 
H = Y^[ E ( k ) - %L+«fc,+ - l E ( k ) + h]a{_a k! _, (23) 
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where E(k) = \JJJ[ - J 2 ) 2 + 4J X J 2 cos 2 k. Clearly, the 
dimerization induces a gap in the dispersion. Using the 
results of the previous sections and Appendix |b[ the ther- 
mal Drude weight of this dimerized chain is 



k XY {T,h) = — 



dk 



(E(k) - hf 



cosh' 



+ - 



{E{ k) + h) 2 

( E(k)+h \ 
\ 2T J 



cosh' 



V dk J 



(24) 



For T -C | Jx - J 2 |, and \ J X - J 2 | < \fJ^h the physics 
is similar to that of the continuum model of the weakly 
coupled ladder discussed in Sec. II C . The fact that the 



model is defined on a lattice allows us to verify that for 
h = and at very high temperatures the thermal conduc- 
tivity indeed decays as T~ 2 . For T ^> V JiJ 2 , and h = 0, 
since the energy spectrum of the XY chain is bounded, 
one has: 

^MT) = ^f"jkE { kf(^-\ (25) 



We note that the integral is independent of the magnetic 
field, so we have a kind of "sum rule" . Since in the pres- 
ence of the magnetic field, it is easily seen that the low 
temperature thermal weight is enhanced by a factor e h / T , 
this necessarily implies that for higher temperatures, the 
thermal weight must decrease when a magnetic field is 
applied This ; " ™„fl,.™^ k„ t?;„ g wt„ 




FIG. 2: The field dependence of the thermal Drude weight 
k(T) for | Ji - J 2 | = 1, and yflUl = 1. 



This result is in fact more general. Since, at high tem- 
peratures, (J 2 ) is finite for a lattice model, we have the 
asymptotic behavior k ~ (J 2 )t=oo/T 2 . Another inter- 
esting limit is when V Ji J 2 <C \Ji — J 2 |, i.e. when the 
spin gap is much larger than the bandwidth of magnetic 
excitations. In this case, for v 'Ji J 2 <C T <C | Ji — J%\, 
replacing E(k) in (|5|) by | J\ — J 2 | = A we obtain 



k{T) 



aT 2 cosh 2 (A) 



(26) 



Note that the Drude weight can be recast in the form 
k(T) = irC v (T)v e g., where C V (T) is the specific heat of a 
fermion that can occupy two levels separated by | J\ — J 2 1 
and with an effective velocity w e //. ~ J\Jij\J\ — ^l- 
However, this analogy cannot be extended systematically 
to other spin gapped systems. 

Turning to the effect of the magnetic fie ld, the Drude 
weight can again be rewritten as (see Eq. (A9)) : 



note that for high magnetic fields, a double peak struc- 

tUre appearS IP ^1t^ ftior-mal ijn^irrVi^ QC QOOn m Tn^irr T'tio 




FIG. 3: The double peak structure in k for h = 0.8, |Ji — Ja| = 
1 and yJuh = 1. 



k{T, h) 



dk[C v {e{k) -h) + C v (e(k) + h)]v 2 {k){27) 



where C v (e) is the specific heat of a single fermion of 
energy e and the velocity v(k) = ede/dk. This form helps 
us derive a kind of sum rule for the thermal conductivity. 
For a free fermion, one has: 



C V (T) 



T 



dT = S(T = oo) - S(T = 0) = k B In 2, (28) 



And thus: 



k(T, h) 
T 



dT = k B In 2 J dkv 2 {k). 



(29) 



double peak results from the low temperature shift of the 
maximum of the contribution of the up spins, in a region 
in which the contribution of the down spins in negligible, 
and the high temperature shift of the maximum of the 
contribution of down spins. A similar double peak is also 
visible in the heat capacity. It would be interesting to 
investigate whether such a double peak is also present in 
other spin gap systems. It is known that a double peak 
is present in the. specific heat of zig-zag spin ladder in a 
magnetic field.Ej To summarize, the sum rule (|29| ) for the 
thermal Drude weight, holds for all spin systems which 
can be described by an effective theory of non-interacting 
fermions. 
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E. The Massive Boson Model 

We now consider the massive triplet boson model 
which was proposed as a phcnomcnological model for the 
spin-1 Heisenberg chain. This ..model can be obtained 
from the nonlinear sigma modelcJ that describes integer 
spin— S chains in the limit S — > oo by softening the con- 
straint on the 0(3) fields. More precisely, this model is 
characterized by a Hamiltonian densitycil 



(30) 



a=l 

where V(<j>) = ^cf> 2 + |(0 2 ) 2 , and [<f> a (x), U^x')] = 
iS(x—x')5 aj ff. The energy current takes the simple foruO 



3e(x) = 



-U 2 ^2 Rad x (f) 



(31) 



Note that this current is independent of the potential 
V(4>) and up to a prefactar, it is just the momentum 
density of the boson fieldEI Consequently, translation 
invariance implies that the total thermal current J e is 
conserved. This allows us to use the Eq. (||) to obtain 
the thermal Drude weight. Since the bosons are weakly 
interacting, we can consider the case A = 0, to obtain the 
Drude weight 



k(T) 



3 
8T 5 



41,2 



dk- 



u k 



sinh" 



2 / ^Qfc) 2 +A 2 



IT 



(32) 



As before, the limit T <C A again leads to the result 
(O) for k(T) and for T > A, we recover a linear weight 
k{T) = n 2 uT. We note that in the case of the Takhtajan- 
Babujian spin-S* chains (whick-are described by SU(2)2S 
WZNW models at low energyc3) , the weight is given by: 



R(T) 



35' 



3SH 



uT, 



(33) 



for S — > oo, this weight is the same as the one of the 
triplet of bosons. This is consistent with the fact that 
the non- linear sigma model describes spin— S chains in 
the limit S —> oo. 



F. Discussion 

In the preceding sections, we have seen that all the 
spin gap systems studied in this paper exhibit the same 
generic behavior for the thermal conductivity. The rea- 
son is that for gapped ID systems that can be bosonized, 
the low energy theory is Lorenz invariant, and excitations 
are described by massive particles having relativistic-like 
dispersions e a (p) — yj (vp) 2 + m 2 a , with the gap A being 
the mass of the lightest particle. When these excitations 



are spin triplets, the lowest excited state contains ex- 
actly one of these particles, and the total energy current 
is e(p) J| which then yields a weight: 



R{T) 



/**>■ (I 



(34) 



Since Lorentz invariance dictates that e(p)§| = P, one 
obtains the same thermal weight as in (|l8|) in the low- 
temperature regime. Examples of systems possessing; 
this triplet branch are the_alternating spin-i chainBc3, 
the two-leg spin laddeiMc3 and the Heisenberg spin-1 
chained. We therefore, expect that the above mentioned 
systems will exhibit a finite thermal Drude weight. How- 
ever, this result could differ in the case of the zig-zag 
ladder or the frustrated spin | chain. This stems chiefly 
from the fact that though the zig-zag ladder has a gapful 
spectrum, the low energyj-espitations having a relativis- 
tic dispersion, are spinonsEIEj carrying a spin | . Another 
example with spinonic excitations is the XXZ chain in 
the Ising phasea. Since the total spin of the system can 
only vary by an integer, the spinons occur in pairs. Con- 
sequently, the interaction between these spinons has a 
strong influence on the thermal weight. In the case of 
the XXZ chain, since the spinons are non-interacting, 
the current of a given excited state is conserved, and 
one expects to recover a finite Drude weight. However, 
in the case of the zig-zag ladder or the frustrated spin 
i chain, the interaction between the spinons can lead 
to a non-conservation of the current of the two spinon 
state, Resulting in the suppression of the thermal Drude 
weightO. 

It would be worthwhile to compare our predictions 
for the Drude .weight for various systems with numeri- 
cal simulationsllZI or with other analytical techniques on 
the lines of Ref. ||l] in the case of integrable models. How- 
ever, in the former case, the extraction of the power law 
prefactor in the activated thermal Drude weight from nu- 
merical data might prove very difficult. 



III. EFFECT OF IMPURITIES 

We have seen in Sec. || that the thermal conductiv- 
ity in clean systems has a Drude peak as a result of the 
translational invariance of the system. In a real system, 
we expect this Drude peak to be replaced by a finite 
thermal conductivity, due to the finite lifetime of eigen- 
states of the Hamiltonian induced by phonon or impurity 
scattering. In the present section, we study the effect of 
impurity s catte ring on the gapped systems we discussed 
in sections II C and [I D . We will begin with a calculation 
of the conductance of the system with a single impurity, 
and then we will turn to a system with a nonzero con- 
centration of impurities. 
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A. Single-impurity problem 

The thermal conductivity of a system with a single- 
impurity xan be calculated using the simple Landauer 
approachES provided, the elementary excitations are non 
interacting. The basic ideauS is to consider two reservoirs 
at temperature Ti and T 2 (with Tj > T 2 ) i n presence of 
a barrier (the impurity potential). Reservoir 1 emits a 
particle with momentum k > 0, energy e(k) and velocity 
de(k) I dk. The probability to traverse the barrier is given 
by the square of the transmission coefficient |i(fc)| 2 . The 
current flowing from reservoir 1 to reservoir 2 is: 



2ir 



Ok 



(35) 



and similarly the current flowing from to reservoir 2 to 1 



be fermionized, so that the resulting low energy Hamil- 
tonian of the ladder with a rung defect reads: 



H 



+i 



0=0*' 

J dxm(x) (X)SS-3SeS 



(41) 



where m(x) = m + g8(x), with to = Jj_/(2n) and 
g = (J'j_ — J±)a/(2n). Clearly each Majorana mode is 
scattered independently from the barrier, so that their 
contributions to the thermal conductivity is additive. As 
discussed in in the preceding sections and in App. |X[ we 
use the correspondence between the Majorana and Dirac 
fermions to calculate the thermal conductivity with the 
barrier. The First Quantized Hamiltonian for the Dirac 
fermions reads: 



J* 



^n 2 (k,T 2 )\t(k)fe(k) delL) 



2tt 



Ok 



(36) 



where ^1,2(^)^1,2) are the fermion distribution functions 
at temperature Ti,2- In the limit T\ ~ T 2 , the net current 
flowing through the barrier is: 



J — ^l->2 — Jl^l — 



/C(Ti)(Ti-T 2 ) 



dk \t(k)\ 2 e 2 (k) Ti - T 2 de(k) 



27r 4cosh (£«; 



T 2 



dk 
(37) 



Hence a knowledge of the transmission probability \t 2 \ 
permits us to obtain the thermal conductance 



K,(T) 



e 2 (k) dejk) 
27r rvvi 4T2 C osh(^) dk ' 



^\t(k)\ 2 



(38) 



We now apply this general formula ( pq) to two spin gap 
systems in which the elementary excitations are non- 
interacting. 



1. Ladder with a defect 

We consider a two leg spin 1 /2 ladder with a defect on 
a rung, described by the Hamiltonian: 

H = J|| Sn tP -S n +l,p + J_L S W| 1 -S„ 2 + Jj_Sq ; 1-Sq ; 2 

(39) 

This Hamiltonian can be fermionized following Ref. pU 
The perturbation to the ladder becomes: 



(J' ± - J 1 _)a 2 (3 1 + m)(0) • (J 2 + n 2 )(0), 



(40) 



where Ji 2 and n 12 are the uniform and staggered spin 
densities, respectively, and the most relevant contribu- 
tion is (J' ± — Ji)a 2 ni(0) • n 2 (0). This contribution can 



H = —iva^dx + m(x)a 2 , 



(42) 



where <Xj are Pauli matrices. Solving the Schrodinger 
equation with appropriate boundary conditions for the 
wavefunction at the barrier, we obtain the transmission 
probability 



\t(k)\ 2 = cos 2 tb 



k 2 + K 2 ' 



(43) 



where K — m/vsinip and tan(^) = 
J_l)/(27t 2 J||), Eq.@ becomes: 

k 2 

(vk) 



JL 

2v 



= (A 



f°° dk 

JC(T) = / f- cos 2 ^ 
Jo 27r 



{vk) 2 



fc2 + if2 4T 2 cosh(^i±^) 
3, (44) 



where we have used e(k) = yj (vk) 2 + in 2 . In the limit 
T — > 0, the transmission probability is dominated by 



momentum k <C K for which |<(fc)| 2 
barrier is a strong scatterer, and 



k 2 /K 2 



1. e. 



the 



K(T) 



3m 
2^T 



e- m / T cot 2 V, 



(45) 



where we have taken into account the triplet of Majorana 
modes . One can obtain an estimate of the temperature 
T* below which this result is valid by noting that for 
T <C to, one has ((vk) 2 ) = mT, so that the criterion 
for low temperature is mT -C to 2 sin 2 ip, i. e. T <C 
T* = J±(J' ± - J±) 2 /J 2 . This temperature is clearly 
much smaller than the gap m. For higher temperatures, 
T* -C T <C to, the thermal conductance is obtained by 
making the approximation |t(fc)| 2 ~ cos 2- 0, leading to: 



AC(T) 



3 cos 2 ib 

cos 2 ib 
+ — — - 

47T 



T I dx 

m/T 



x/2 



cosh(x/2) 



T 



dx 



x/2 



3m/T V COsh (^/ 2 ). 



(46) 
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For T* < T < m, one finds JC{T) ~ tfye- m / T , and 
for T > m, JC(T) = ^f- cos 2 ip. Contrary to the result 
( |l8| ) for the pure ladder, the thermal conductance ( f45| ) 
for T <C T* is purely activated without any T dependent 
prefactor. Therefore, the Drude weight in the thermal 
conductivity for the pure system is not an accurate indi- 
cation on the behavior of the thermal conductivity in a 
system with impurities. The reason for that is clear from 
( p3| ) , namely low energy modes experience much stronger 
impurity scattering than the high-energy ones. It is only 
in the high temperature limit T 3> m that the replace- 
ment S(u>) — ► t is justified. We will see in the following 
section that this result is not restricted to the spin ladder. 



2. XY chain with a defect 

We consider again the XY-chain with alternating ex- 
change of Sec. II D| . We now suppose that the bond 
strength Ji between the sites and 1 is replaced by J[. 
This bond acts as a barrier and using the results of Ap- 
pendix the transmission probability across this barrier 
is given by 



\t{k)f 



4J 2 (J() 



(J 2 -(J{) 2 ) 2 +4J 2 (J{) 2 sin 2 k 



(47) 



In particular, we can show that when k ~ 7r/2 we have: 



l*(*)| 2 = 



WJ?(J[) 2 Jl(k- ir/2) 2 



(J 2 - (J() 2 ) 2 (Ji - J 2 ) 2 + 16J 2 (J{) 2 J|(fc 



(48) 



which indicates that for low temperatures T <C yJx J 2 , 
the behavior of the thermal conductivity in the XY chain 
with a bond defect is identical to the behavior of the ther- 



mal conductivity in the ladder discussed in Sec. Ill A 1 



For high temperatures, T 3> V J1J2, we can neglect the 
variation of the transmission coefficient with the energy, 
and assume that all states have the same probability of 
occupation. Then, the thermal conductance reads: 

AJx+Jz)/T 

JC(T) = T / dee 2 {\t{e)\ 2 ) ~ 1/T 2 (49) 



B. Many impurities case 

In this section, we consider the effect of a finite concen- 
tration of impurities on the thermal conductivity of the 
ladder. As before, the disorder we consider is a random 
rung coupling. The Hamiltonian of the disordered ladder 
reads: 



H 



^„E 



(50) 



where J]_ = Jx + Vi- We have \r/i\ < J±, so that all 
rung interactions remain antiferromagnetic. This Hamil- 
tonian can be analyzed h\ mapping onto a random mass 
Majorana fermions modelEEl. 



H = - 



4 

E 

a=l 



-IV , 



dx{-[i a R {x)d x e R {x) - Cdx)d x C L {x) 



+im a (x)e R (x)a(x)}}, (51) 

with m 1 ' 2,3 (x) = m(x) for the triplet magnetic exci- 
tation, m° = — 3|m|(a;) for the singlet excitation and 
m(x) = m + r/(x) where rj(x)rj(x') — D5(x — x'). We 
note that disorder does not mix the different flavors of 
Majorana fermions. Consequently, the contribution of 
the Majorana modes to the thermal conductivity remains 
additive. As before, to calculate the disorder induced 
self-energy, it is useful to recast the above problem in 
terms of Dirac Fermions. 



H = -iv I dx(ij R d x ip R - ip' L d x ip L ) 



+m(x) / dx(ip R ipL + 4> l iPr), 



(52) 



We note that the Hamiltonian (^) can also be, derived 
from a dimerized XY chain with bond defects. E 2 ] We de- 
fine the 2x2 matrix disordered averaged Green's function 
G by its components, 

= -<T r ^ a (a:,r)^J(0 ) 0)) J (53) 

where a, (3 £ {R,L}. The Hamiltonian ( |52] ) can be 
rewritten in matrix form as: 



H 



dx^ (x)[-ivT 3 d x +m(x)n]^!(x), (54) 



where rx t 3 are Pauli matrices. The impurity self-energy 
matrix S can be calculated -within the Self Consistent 
Born Approximation (SCBA)EJ and satisfies the Dyson 
equation for the disorder averaged Green's function: 



(iujn — vkr^ — mn — S)G = 1, 



(55) 



In this approximation, the self-energy is independent of 
momentum and is determined self-consistently by 



f dk 

£(mj„J = D / — Ti[iuj n - vkr 3 - mn 

J 27T 



(56) 

Clearly, the self energy possesses the following structure, 
T,(iu! n ) = ia(iu> n ) + V(iu> n )Ti leadipg to the following 
self-consistent equations for a and VE3: 



a(iuj n ) = 
V(iw n ) = 



D 



a(iuj n ) - LU n 



2v ^{(j(iuj n ) - uj n ) 2 + (m + V(iuj n )) 2 
D V(iuj n ) — m 



(57) 



2v y/((r(ibj n ) - oj n ) 2 + (m + V(iuj n )) 2 

Introducing the dimensionless variables: s = ia/m, 
t = V/m, x = ibj n /m and A = D/(2vm) the above self- 
consistent equations simplify to: 

(58) 



2s' 



s 3 x-s 2 ( 



1 



A 2 

A 2 ) - X 2 sx+ — x 2 = 0(59) 
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ReP(e + iO,e + iO) in Eq. (60) can be neglected, we find 
that this term is indeed crucial to take into account the 
presence of a gap in the energy spectrum. To proceed 
with the calculation of k, we first note that for weak dis- 
order i.e., D <C vm, since the off-diagonal self-energy V 
always occurs in the combination m + V (|55|), it is rea- 
sonable to neglect V in the Green's function G which can 
then be approximated as 



G(k,Lu) = 



LO — u(lj) + vk<73 + 771(7l 

[lj - a(to)} 2 - e(k) 2 ' 



(62) 



FIG. 4: The imaginary part of the self energy in units of 
m calculated from Eq. (B7]) for various values of the scaled 
impurity strength A. Note that the effective gap decreases 
with increasing A. 



This quartic equation can be solved numerically. Some 
sample curves are shown in FigJ^. We find that the dis- 
order renormalizes the gap in the spectrum and a suffi- 
ciently strong disorder (A = 0.5) closes the gap indicating 
a disorder induced phase transition within the SCBA. A 
plot of the renormalized gap lj c as a function of disorder 
strength is shown on Fig. |5| In the ensuing calculation, 
we only consider disorder strengths for which the renor- 
malized gap is 




where e(k) = y/v 2 k 2 + m? . This yields 
1 vka3 + mo\ s 



G(k,u + i0) 



1 



G(k,uj-i0) 



2 2e(k) ) LO- u(lj) - e(k) 

f 1 vkoz + mo\ \ 1 

^2 2e(k) J cj-a(u) + e{ky 

1 vkcTz + mo\ \ 1 

2 + 2e(fc) ) uj-a*(Lj)-e(k) 
( 1 vkus + mo\ \ 1 

+ V2 Mk) J uj-a*{uj) + e(k)' 

Substituting the above in (|6l"|), we obtain 

dk 

(63) 



P{u}+iO,u-iO)-ReP(uj+iO,LU+iO) = I ^{v 2 k) 2 K{u,k) 

ZTT 



where 
K(w,k) 



(Imcr(w)) 5 



[(lj - e(k) - R C( t(cj)) 2 + {lmcr(Lu)) 2 } 

(Imcr(w)) 2 
[(w + e(k) - Recr^)) 2 + (Imcr(w)) 2 ] 



r(64) 



/ rrt\ 1 f duj ( dn F 



where: 



This expression can now be used in ( J60| ) to obtain the 
thermal conductivity. At low temperatures, the deriva- 
tive of the Fermi function in ( p0[ ) decays exponentially as 
cxp — (u>/T) indicating that frequencies much larger than 
T can be neglected in the integral for the thermal conduc- 
tivity k(T). Consequently, the low temperature behavior 
of k(T) is completely dictated by the the low frequency 
behavior of K (u>, k). We now analyze the behavior of K. 
Firstly, since the diagonal self energy Imcr(w + iO) = 0, 
for u) < lo c i.e., for frequencies smaller than the disor- 
der renormalized gap, K(u>, k) is identically zero for all 
lj < lo c . A typical plot of if as a function of to for two 

[P(Lj-iO,Lj+iO)-RcP{oj+iO, cj+i^ eTent values of k is shown on fi S ure I- 

Clearly, the dominant contribution to k for temperatures 
T < lj c comes from the behavior of K (lj, k) in the vicinity 
of lo c . This behavior has been analyzed numerically, and 
we find that for all k, K can be developed as a series in 
lj — Lj r : 



FIG. 5: The effective gap in units of m as a function of the 
scaled impurity strength A. 



Using the results of the previous sections, the thermal 
cogluctivity for the disordered ladder can be rewritten 



(60) 



/dk 
|^(7; 2 fc) 2 Tr[G(fc,c)G(fc,^)], (61) 

Vertex corrections to ( |60| ) are negligible in the weak dis- 
order limit. Contrary to the suggestion in Ref. Efl that 



K(uj,k) = a(k)(uj — lj c ) 

+P(k)(uj — to c ) 2 + 



(65) 
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K(to.k) 40 




FIG. 6: The function K(u>, k) in units of m 2 for A = 0.1 and 
various values of k (in units of m/v.) 



We find that a, (3 are fast decreasing functions of |fc|, such 
that J dkk 2 a(k) < oo and J dkk 2 (3{k) < oo. Substituting 
(p5|) in (|60|), we obtain 



K (T) = J dk J duj fc 2 e~" /T e(w - w c )[o(Jfc)(w - lj c ) 

+P(k)(LJ-LJ c ) 2 + ...] (66) 

leading to, 

k(T) ~ &e-^ /T + fiTe-^^ + o(Te~^ /T ) (67) 

a, (3 are temperature independent constants. We find 
that the results are similar to those for a single impurity 
( [l5| ) with the difference that a finite concentration of im- 
purities renormalizes the gap. At very high temperatures 
one recovers the usual T~ 2 decrease of the thermal con- 
ductivity. It would be interesting to obtain the cross-over 
behavior from the the low to high temperature regimes. 
However, the rather complicated form of the self energies 
makes analytical calculations very difficult for these in- 
termediate temperatures and this is left for future work. 
For small magnetic fields h <C ui c the same result holds 
with the substitution lu c — ► uj c — h. 



We see that in the limit of very small self energies, we 
recover the Boltzmann equation result for the thermal 
conductivity!^. Comparing ( |6^ ) with (||) and (17), we 
see that if Ima(e(k)) = t~ x where r is a constant in- 
dependent of e, the thermal conductivity can be writ- 
ten as a product of the Drude weight in the absence of 
impurities and the mean relaxation time r as was pro- 
posed in Ref.|l5|. The underlying assumption there, was 
that all the eigenstates of the system have the same life- 
time r independent of the energy of the eigenstate. As 
shown above, the explicit energy dependence of the self 
energy found in Sec. [II, i.e., the fact that the low energy 
spin excitations are much more scattered by impurities 
than high energy excitations, shows that any assump- 
tion of energy independent lifetime is invalid even for the 
simplest models. As a result, the thermal Drude weight 
can atmost yield a heuristic behavior of the thermal con- 
ductance in a real system in which spin excitations are 
interacting with impurities and/or phonons due to the 
different extrinsic lifetimes of current carrying states. 

We present a brief comparison of our results with ex- 
periments. Since disorder is ubiquitous in real systems, 
it is reasonable to compare our results for the two leg 
ladder with impurities with experimental measurements. 
One of the systems studied extenspflljjja the spin gap 
compound Sri4_ x (La, Ca) x Cu2404ili3lI3ii-J. In the insu- 
lating phase, these systemSr-can be well described by an 
array of two-leg spin ladders^. In this system, the phonon 
subtracted thermal conductivity was shown to have an 
exclusive spin contribution. At low temperatures, a fit 
for spin thermal conductivity yielded a k(T) ~ e~ A / T . 
This low temperature fit is in good accord with our pre- 
diction of k(T) ~ e~ Wc l T . However, a full quantitative 
comparison requires an understanding of the effect of the 
disorder in the material on the spins, the effect of spin- 
phonon interactions and other exchange interactions in 
the ladder. 



CONCLUSION 



IV. DISCUSSION 

We now highlight the connection between our approach 
and that of the Boltzmann equation. A Boltzmann like 
equation for the thermal conductivity can be recovered 
from the SCBA in the limit of small self energies. For 
a(ui) <C m, the function K defined by fl64|), takes the 
simpler form 

K(uj, k) = tt[S(uj - e(k)) + 5(uj - e(fc))](Imcj(w))" 1 (68) 
Inserting the above result in (|60|), we obtain: 

(69) 



In the present paper, we have calculated the thermal 
conductivity of gapless spin chains and spin gap systems 
including the two- leg spin-i ladder and the dimerized 
XY spin-i chain. In the absence of disorder, the ther- 
mal Drude weight of gapless spin chains vanishes linearly 
with temperature. On the other hand, for the ladder 
and the XY chain, which can both be represented as 
free fermions, we find that the thermal Drude weight 
k ~ j>— i/2 e -m/T^ -y^hgj-g m ^ the gap to the lowest 

triplet excitation. For intermediate temperatures, k«T 
and decays as T~ 2 at very high temperatures. We ar- 
gue that this behaviour is generic to all quasi-one dimen- 
sional spin gapped systems having low energy triplet ex- 
citations with a relativistic-like dispersion. We have also 
considered the effect of a magnetic field which results in 
a substantial enhancement of the low temperature ther- 
mal conductivity. Furthermore, in the case of dimerized 
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XY chains, a double peak can be obtained in the ther- 
mal conductivity for large enough fields. We have also 
studied the effect of impurities on spin gap systems like 
the ladder and the XY chain. Impurities destroy the 
Drude peak resulting in a finite thermal conductivity at 
zero frequency. This thermal conductivity has a generic 
form e~ m / T at low temperatures where, to is the effective 
gap of the system. It would be interesting to include the 
effects of magnetic impurities and also scattering from 
phonons. These and other questions are left for future 
work. 
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us to rewrite the total energy current J e = J dxj e (x) = 
Efc v2 Kcl,R c k,R. + c l,L c k,L ~ 1) as: 
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Using (A4), one easily obtains 



(J e 2 ) = 2]>> 2 fc)> + (fc))(l - (n+(k))) (A6) 



where the Fermi distribution function (n + (k)) = (e' 3 ' 5 ^)- 
From (||), the thermal conductivity k(u,T) - 
k6(lu) with a Drude weight 
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APPENDIX A: CALCULATION OF THE 
THERMAL CONDUCTIVITY FOR MAJORANA 
FERMIONS 

1. Thermal conductivity of Dirac fermions 

We consider massive Dirac fermions with the following 
Hamiltonian density: 



1 f c 



41,2 



dk- 



v k 



cosh' 



(#) 



(A7) 



The above result has a simple interpretation in terms of 
kinetic theory. Consider the expression for the specific 
heat: 



C V (T) 



1 

J2 



dk 



<kf 



2tt 



cosh' 



(#) 



dk 
2^ 



c v (k,T), (A8) 



i.e., a mode of momentum k contributes c v (k) to the spe- 
cific heat. Such a mode has a velocity v(k) = v 2 k/e(k). 
This now permits us to rewrite (16) as: 



H(x) = -i-{i^ R d x ilj R 



+m{e^i, R il> L + e-^i>li }R ), 
In this case, the energy current reads: 



(Al) 



dk 
2^ 



c„(/c,T)u 2 (fc), 



(A9) 



i.e., the contribution of each mode k to the the Drude 
weight is just the product of its specific heat and square 
of the velocity. This is similar to the kinetic theory result 
i~{dx4> R )i>R.+^ L d x 'il>L-{dx^ L )'>pL)(x) that the thermal Drude weight is given by the product 



j e (x) = -i—{tp R d x ^ 

We note that using the transformation ipi — > fipLi we 
can reduce the Hamiltonian flAip to the case = 0, 
while leaving the current ( |A2|) invariant. Therefore, for 
the purpose of calculating the thermal transport we can 
without loss of generality restrict to the case <f) = in 

Using the Fourier decomposition 



of the specific heat and the square of the velocity of the 
free modes. 



Majorana fermions 



(A3) 



It is well known that the Dirac Hamiltonian (Al) can 
be re-expressed in terms of two Majorana fermions fields 
defined by ip v — {Q, + iQ,)/\f2. The Hamiltonian can be 
written as a sum of two Majorana Hamiltonians 



the Dirac Hamiltonian H = J dx7i(x) can be diagonal- 
ized to obtain 



Hduhc — HmIC 1 } + Hm[( 2 ] 



(A10) 



(A4) 



where the fermionic operators Ck.± are linear combi- 
nations of the c k R j L such that c£ + Ck,+ + c|, _Cfe,_ = 
c\ R Ck,R + c\ L Ck,L and e(k) = yj (vk) 2 + to 2 . This allows 



Similarly, the energy current (A2) can be written as the 
sum of two energy currents, each associated with one Ma- 
jorana field: j®(x) = jl(x) + j 2 (x). The thermal con- 
ductivity of the Dirac Hamiltonian can then be written 
as the sum of the conductivities associated with the two 
Majorana field i.e., KDi rac (^,T) = k 1 (w,T) + k 2 (uj,T). 
The expression of the currents and the Hamiltonian be- 
ing identical for the two species of Majorana fermions, it 
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is clear that k 1 (w,T) = k 2 (uj,T). Thus, one obtains the 
generic result that 



KMajorana(w, T) — KDirac (w, T)/2. 



(All) 



This result shows that it suffices to calculate the thermal 
conductivity of the corresponding Dirac Hamiltonian to 
obtain the Majorana thermal conductivity. This corre- 
spondence holds provided there are no interactions be- 
tween the various species of Majorana fermions. 



APPENDIX B: THERMAL CURRENT IN THE 
PRESENCE OF AN APPLIED MAGNETIC FIELD 

In the presence of an applied magnetic field, the Hamil- 
tonian density is TC(x) = Ti.(x) h=0 — h ■ m(x), where 
m(x) is the magnetization density and h is the exter- 
nal magnetic field. Using the continuity equation for the 
Hamiltonian density and the equation of conservation of 
the moment ftm + d x j s = 0, the thermal current is now 
given by 



jth.( x ) = 3e{x) - h • j s (x), 



(Bl) 



where j e is the energy current for h = and j s (x) is the 
magnetization current. 

For the specific case of the ladder with a magnetic field 
along the z direction, the Pauli coupling is 



fM&& + &&) (B2) 



Note that the contribution to the thermal conductivity 
arising from the £L L , L is not changed by the applica- 
tion of the magnetic field. To obtain the thermal conduc- 
tivity coming from the, modes £ 1,2 it is convenient to turn 
to the Dirac FermionsBE3 ip u>a = (S + *£5)/V2. Then, 
one can rewrite H ma „ as: 



h / dx(ipt s ip R<s + Vx, a i/>L, a ) 



(B3) 



The expression of the total thermal current then reads: 
T de 



(B4) 



The contribution of the £ 1,2 modes to the Drude weight 
in the thermal conductivity is then calculated to be 



k (T, h) + k\T, h) 



1 

8T 2 



dk 



cosh ^W-^ 
(e(k) + h) 2 



cosh 



(B5) 



APPENDIX C: EIGENVALUES AND 
EIGENSTATES OF THE FERMIONIZED XY 
CHAIN 

1. Translational Invariant Case 



The eigenstates of the Hamiltonian (E2J) are obtained 
by solving the equations: 



J\Al n + J 2 ^4 2n +2 — EA 2n +i 
JlA<ln—\ + ^1^2)1+1 = EA 2n 

One finds positive energy solutions: 



A 2n 

A 2 n+1 



2 lkn I e" #fc/2 
~~ ' 1 J>i>k/2 



(CI) 
(C2) 



(C3) 



with E(k) = y/(,h - J 2 ) 2 + 4Ji J 2 cos 2 k and J x 
J 2 e 2jfc = E(k)e lc ^ k , and negative energy solutions: 



A 2n 

A 2 n+1 



2ikn 



e -i<Pk/2 



(C4) 



with E(k) = - 1/( Ji - J 2 ) 2 + 4J\ <h cos 2 k and J\ + 
J 2 e 2ik = \E(k)\e i,l>k , 

2. Single impurity case 

Clearly, the solutions with momentum k and — k are 
degenerate in energy, thus we search the solution as a lin- 
ear combination of these solutions. The system of equa- 
tions to solve reads: 



J 1 A 2n + J 2 A 2n+2 = EA 2n+ i(n 0) 
J 2 ^2n-i + JiA 2n +i = EA 2n (n j£ 0) 
J[Ao + J 2 A 2 = EA 1 (n = 0) 
J 2 A^ + j[Ai = EAo(n = 0) 

We search for solutions of the form: 



(C5) 
(C6) 
(C7) 
(C8) 



A2n \ 2/7. n 

A2n+l 



-i<j> k /2 
i<Pk/2 



for n < — 1 and: 
A 2n 

A 2 n+1 



= t(k)e 



2/7/. 



-r(k)e 



-i<f> k /2 
i<Pk/2 



— 2ikr< 



e i<f>k/2 
g— #fe/2 

(C9) 



(CIO) 



for n > 1. Applying equations ( jCq ) for n = — 1 and ( |C6| ) 
for n = 1 we obtain respectively: 



A Q = e -«t>k/2 + r ( k y4>k/2 

A x = <(/c)e i0fe/2 



(Cll) 



The equ ations that determine t, r are obtained from (C7) 
and (C£). They read: 



J 2 ( e - 2l V^/ 2 + r(fc)e 2j ' £ e~^ fc/2 + J[t{k)e** k ^ = £;(e- 1 ^/ 2 + r(i 
jJ( e -***/a + r (fc) e ^/ 2 ) + J 2 t(k)e 2lk e- l ^l 2 = £t(fc)e^ fc/2 
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Using the relation Ji + J 2 e l2k = E(k)e l<t ' k these equations 
are simplified as follows: 

- Jie-^ /2 r(fc) + J{e^ fc/2 t(fc) = J 1 e l,pk/2 (C14) 
J' 1 e l4 ' k/2 r{k) - J ie -^/ 2 t(fc) = - J[e- l<t »</ 2 (C15) 

We obtain the transmission amplitude t(k) and the 
reflection amplitude r(k) as: 



t(k) = 



—2iJ 1 J[ sin 4>k 



(C17) 



r(k) = ^ ~ 31 
[ ) J 2 e~^ - (J{)V<^ 



(C16) 
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